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Abstract. The leading order coefficients of the /^-function of QCD are computed in a large Ay- 
expansion. They are in agreement with the three loop MS calculation. The method involves 
computing the anomalous dimension of the operator (Cr^,^)^ at the d-dimensional fixed point 
in the non-abelian Thirring model to which QCD is equivalent in this limit. The effect the 
0{1/Nf) corrections have on the location of the infrared stable fixed point for a range of Nf is 
also examined. 
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The strong force of the standard model is described by quantum chromodynamics, (QCD), 
which is an SU{3) gauge theory which is asymptoticahy free. At large energies the fundamental 
quarks behave as though they were non-interacting. In terms of the field theory itself this 
property is a consequence of the leading coefficient of the /3-function being negative. The initial 
one loop calculation was carried out in Higher order corrections have also been determined. 
The remaining scheme independent term, ie the two loop contribution, was calculated in [^. 
The three loop term was computed in the Feynman gauge in S using dimensional regularization 
in the MS scheme. More recently this result was checked in Q in an arbitrary covariant gauge 
where the remaining renormalization group functions, such as the gluon anomalous dimension, 
were also deduced in arbitrary gauge, The three loop quark mass anomalous dimension, 
which is gauge independent, is available too, ||5|, One reason for such precise information is 
that, for example, it allows one to obtain a more accurate insight into the variation of quantities 
with energy scale. Fundamental in this respect is the /5-function as it always appears in the 
appropriate renormalization group equation, (rge). 

These higher order analytic calculations of the rge functions are exceedingly tedious to com- 
pute, however, due to the huge number of Feynman diagrams that arise. For such results to be 
credible it is important to have independent checks on the expressions obtained aside from the 
obvious one of performing another complete evaluation which may be a waste of resources. In 
this letter we provide the results of such a procedure for the QCD /3-function. This is the large 
Nf technique of determining exact all orders results of the rge functions of gauge theories at suc- 
cessive orders in powers of ^/Nf, where is the number of fundamental fields. The technique 
was initially developed for low dimensional models in a series of impressive papers, ^, ^. 
Briefly the method involves computing appropriate critical exponents at the d-dimensional fixed 
point of the /3-function as Nf ^ oo. Through the critical rge these d-dimensional exponents 
encode all orders information on the coefficients of the corresponding rge function. Clearly the 
values will overlap with the lowest known orders, providing the partial check we have indicated. 
From a technical point of view one benefit of this approach is the exploitation of the conformal 
symmetry at the fixed point which simplifies the resummation of the minimal set of Feynman 
diagrams comprising the relevant Schwinger Dyson equation. The calculation of the 0{1/Nf) 
QCD /3-function here completes the leading order analysis as the quark, gluon and ghost dimen- 
sions were deduced in |jlO[ , in the Landau gauge and agreed with the three loop results of ^ . 
Another motivation arises from a comment in Q in regard to future calculations. It is indicated 
that the four loop QCD /3-function is attainable. The main obstacle, though, would appear to 
be correctly generating and treating the vast numbers of Feynman diagrams. Therefore the new 
coefficients we will deduce from our results will be important in this respect. 



We recall that the 0{\/Nf) computation of the QED /3-function is available, |12]. That 
calculation was carried out by inserting the implicit bubble sum of the photon propagator 
in the 2 and 3 point functions and then deducing the MS coefficients of the renormalization 
constants using dimensional regularization. Those results have been reproduced in the critical 



point approach, [13|. One interesting aspect of |12| was the search for other fixed points in the 



strictly four dimensional QED /3-function, for a range of values of the coupling. Although none 
were observed it would be a worthwhile exercise to repeat that analysis in the non-abelian case 
especially as at two loops such a point exists, [Q, for a range of Nf. 

We recall the fundamental ingredients for treating QCD in large Nf in our approach. The 
lagrangian is 

L = - (1) 

where ip'^^ is the quark field, A'^ is the gluon field, the covariant derivative is = 9^ + iT"'A'^ 
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with the generators of the colour gauge group and 



IS 



the field strength with / the structure constants. The ranges of the indices are 1 < i < A^, 1 
<a< {N^ - 1) and 1 < I < iVc and the Casimirs are tr(r'^r^) = T{R)6''\ T'^T'' = C2{R) and 
facdfbcd ^ C2{G)6''\ To three loops, in d-dimensions, [1-4], 
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(d - 4)5 + 
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-T{R)Nf - -C2{G) 



a 
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]^C2{R)T{R)Nf + ^C2{G)T{R)Nf 
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Cl{G) 
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G2{R)T\R)Nf + ^C2{G)T\R)Nf + ^Cl{R)T{R)Nf 
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C2{R)G2{G)nR)Nf 
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■Ci{G)T{R)Nf + 



1728 



<7^ + 0(/) (2) 



where our coupling g \s g = (e/27r)^. The presence of the 0{g) term of (2), corresponding to the 
dimension of the coupling in d-dimensions, gives rise to our non-trivial fixed point, gc- Explicitly 



3e 



9c 



+ 
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T{R)Nf 4T^{R)N^ 
+ (^G2iR) + ^C2{G)]e^ + 0{e^^ 



33C2(G)e - (27C72(/?) + 45C72(G)) 




(3) 



where d = 4 — 2e. In the neighbourhood of this point the quark and gluon anomalous dimension 
were deduced in the Landau gauge as, at leading order in 1/Nf, respectively, [p^], 



ri + X 



C2iR)v? 
T{R)Nf 

G2{GH 
2{pi-2)T{R)Nf 



(4) 
(5) 



where r/° = (2/i-l)(/i-2)r(2/i)/[4r2(/i)r(^+l)r(2-/i)] and d = 2/i. Moreover, the asymptotic 
scaling forms of the respective propagators are, as k"^ — > oo, 



A^t 



fi—a 



A^u{k) 



B 



A;2 



(6) 



where a = /i + ^rj, (3 = 1 — rf — x and A and B are amplitudes though only the combination 
z = A^B appears in calculations. Specifically z = V{^ + l)r?°/[2(2;U - l)(/x - 2)T{R)Nf]. 

One feature which simplifies the fixed point analysis both for computing the 0{1/Nf) cor- 
rections to (5{g) and (4,5) arises from the universality class in which QCD belongs. For instance, 
it is widely accepted that the 0{N) bosonic a model and the 0{N) (p'^ model are equivalent, or 
in the same universality class, at the d-dimensional fixed point analogous to (3). In other words 
critical exponents calculated in either field theory at this fixed point are the same. So if one 
wished to examine the critical behaviour of the three dimensional 0(3) Heisenberg ferromagnet 
to which both are equivalent then either model can be used. Another widely studied equiva- 
lence is the Yukawa interaction and the Gross Neveu model with the same chiral properties, 
|15| . Likewise the Thirring model and QED are equivalent and have been the subject of recent 
interest, |l^, 17|. For QCD the relevant model is the non-abelian Thirring model, (NATM), 



whose lagrangian is 



L 



l.iJ 



2A 



(7) 
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where is an auxiliary field, which if eliminated leads to a 4-fermi interaction, and A is the 
coupling whose dimension is {d — 2) when compared to the {d — 4) of the QCD coupling. It has 
been demonstrated in |18] that it is equivalent to QCD in the large Nf limit. One feature of that 



work was the correct reproduction of the three and four gluon vertices of non-abelian theories, 
which are absent in (7), by integrating out quark loops in the gluon 3- and 4-point functions. 
For our computation the major simplification is the fact that by calculating with (7) at gc the 
resulting exponents, which are universal, will be equivalent to those computed in (1). Then by 
decoding the exponent using (3), we can deduce information on the perturbative structure of the 
rge functions. A test of this argument will be the correct reproduction of the MS coefficients. 
Importantly though we need only consider graphs which are built out of the single interaction 
of (7). 

We now turn to the details of the calculation. Ordinarily one computes uj = — \l3'{gc) by 
considering corrections to (6) in the Dyson equations, ||7|, |T^. Equivalently one can identify the 
composite operator in the lagrangian whose coupling relates to the ordinary coupling constant, 
1^]. Then the anomalous dimension of that operator is related via a scaling law deduced from 
the lagrangian to (3'{gc)- In QCD the appropriate operator is (G^j^)^ as we use a formulation, 
(1), where the coupling is defined in such a way that the 3-point interaction is ip^^T°"ipA'^. Thus 
the canonical dimensions for the fields essentially satisfy the condition for conformal integration 



or uniqueness, |1£]. From the second term of (1), therefore, we have the scaling relation 

= ri + X + XG (8) 

The quantity XG is the critical exponent corresponding to the renormalization of the pure (com- 
posite) operator (G^^,)^, whilst the gluon dimension arises because of the wave function renor- 
malization of the constituent fields of the composite. Thus computing XG gives w from (8). 

To evaluate XG one substitutes the critical propagators (6) into the relevant 0{1/Nf) set of 
Feynman diagrams and determines the residue of the simple pole in A, [^. This regularizing 
parameter is introduced by the shift /? ^ /? — A. The contributing graphs are illustrated in 
fig. 1 and each is computed in the Landau gauge to avoid mixing with [d^A'^)'^, The first 



two graphs occur in QED and have been computed in |13]. Here we note their respective colour 
group factors are C2{R) and [C2{R) — C2{G)/2\. The third graph arises from the cubic term 
of the operator and it and the final graph are purely non-abelian, each having group factors 
C2{G). The computation was carried out by the application of standard techniques for massless 
integrals including integration by parts and conformal methods. Nevertheless several difficult 
subintegrals lurk within the final graph which were tedious to determine. To verify that we had 
obtained the correct values for them we calculated several graphs of which contained the 
same subintegrals and checked that the total expression we computed agreed with the values 
given in |^]. Useful in this and other respects were the packages Reduce |20| and Form 
We note that the values are respectively, 

/i(/i-l)(2/x-l)7?o (4//^ + /i - 9)r?o _ (4^^ _ 2^2 _ 4^ ^ ^^^o 
+ ' + ' 2(/i + l)(2^- !)(/.- 2) 

(4/x6 - 6//5 + 18/ - 67/i^ + 85/i^ - 38/i + 6)r?? 
4(2// -l)(/. + l)(^- !)(/.- 2) 



(9) 



Further at leading order there are no ghost contributions. One can see this by attempting 
to include them in the formalism we will describe later and then observing that the first ap- 
pearance of any contribution is at 0{1/NJ). This feature was also observed in [^] where the 
non-abelian generalization of the CP{N) a model was studied. Indeed we make use of some of 
the observations of ^ here. 
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The final result is 



oj = in -2) 



^ ^ ^ ; 4^2^ _ ^^^^ _ -^^ 



T(i?)iV/ 
(10) 

A final check on (10) is that it correctly reproduces the 0{1/Nf) terms of the three loop (3- 
function of (2), [1-4]. This amounts to the terms with C2{G) as those with C2{R) have been 
verified for QED in |2|. In three dimensions 

1 10C2(G) 
- " 2 " 2,7.^T{R)Nf ^^^^ 

Prom (10) we can now deduce higher order coefficients which will appear in the MS /3-function, 
by carrying out the e-expansion of (10) and using (3). Thus defining the leading order large 
coefficients by a„, 

oo 

m = /3o5' + E«-+i[^(^)^/]V+' (12) 



n=l 



with /3o = [2T{R)Nf/^ - llC2(G)/6], then 
[154C2(i?) + 53C2(G)] 
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05 



3888 

[(288C(3) + 214)C2(fl) + (480C(3) - 229)C2(G)] 



31104 

«6 = ^^^[(864C(4) - 1056C(3) + 502)^2 (i?) + (1440C(4) - 1264C(3) - 453)^2 (G)] 

«7 = i^^L. [(3456C(5) - 3168C(4) - 2464C(3) + 1206)^2 (i?) 
IbZyblb 

+ (5760C(5) - 3792C(4) - 848C(3) - 885)C2(G)] (13) 



Having obtained the set {an} for non-abelian theories we can examine the purely four di- 
mensional /3-function and search for fixed points other than the well known infrared stable point 
of Banks and Zaks, g^^ ■, [0]. Its existence is important for recent developments in supersym- 
metric theories in relation to electric-magnetic duality, Those infrared fixed points are 

determined by using exact non-perturbative arguments in the limit Nc, Nf ^ 00 with Nf/3Nc 
held fixed, (The one loop coefficient of the /3-function for that model is {Nf — 3Nc) for 

SU{Nc)-) Further in the context of 1/A^ expansions a (16^ — Nj) expansion from g^^ has 
been used to obtain an estimate for as for low Nf, |25 

corrections have on P{g), we first recall properties of g^^ ■ In it was observed that for a 
range of Nf the two terms of the two loop /3-function have a different sign which therefore gives 
rise to a non-zero critical coupling, g^^ ■ For SUc{S) this range is 8 < < 17, |]14|, and we 
have recorded the explicit values of g^^^ in this case in our notation in table 1. Subsequently 
one can also study the effect that the inclusion the three loop term of the /3-function has on the 
location oi g^^ . We have analysed (2) numerically and determined the three loop values of gj^^ 
which are given in the second column of table 1. Several features are apparent. First, the range 
of Ay for the existence of such an infrared fixed point is extended to 5 < Ay < 17. Second the 
effect the three loop correction has is to move the location of g^^ towards the origin. In other 
words to a region where perturbation theory would be valid. These observations, however, ought 
to be qualified. It is not clear whether this picture is meaningful because as Nf decreases g^^ 
clearly increases away from the region where perturbation theory is useful. In other words for 
low values of Ay we can not make a reliable statement on even, say, the three loop range of Ay 
for which g^^ occurs. One indication of where the perturbative picture may not be valid can be 



Before studying the effect 0(1/Ay) 
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deduced from the values of the critical exponent which is a physically meaningful and 

calculable quantity. In table 2 we have given the corresponding values for P'ig^^) as deduced 
from the two and three loop values of g^^ respectively. Clearly the three loop corrections do 
not affect the two loop values appreciably for = 14, 15 and 16 suggesting that higher loop 
corrections are small. For lower values the divergence is evident indicating that the four and 
higher loop contributions would be needed to make an accurate estimate of the exponent. In 
light of the critical coupling being smaller for a larger range of it would better of either, 
though, to take the three loop values of P'{gc^) as being the more reliable. 

Now we consider the effect that the 0{\/Nf) corrections of (10) have. We have studied 
the case A^c = 3 in various ways. First, as in ||l^ we examined the /3-function given by just 
taking all the leading order coefficients a„ which was improved by non-abelianization, [p^. 
This entails replacing Nf by the one loop /3-function coefficient through the shift Nf {Nf — 
llC2(G)/[4r(i?)]). It turns out that in searching for zeroes of the four dimensional /3-function 
that the contributions from these O{l/(3o) coefficients on their own are not sufficient for even 
obtaining a fixed point g^^. This is the same as was found in the QED case, |12], in the range 
of couplings where the series was convergent. Instead, to improve this situation we took the 
two and three loop /3-functions of (2) and then included all subsequent information included in 
(10). The point of view being that one can at least study the effect the 0{1/Nf) corrections 
have on the fixed point g^^ which is known to exist. It turns out that in this approach we 
did not observe any non-trivial fixed points other than g^^ in g > which was independent 
of the number of terms included. From a practical point of view in our analysis we truncated 
the series for f3{g) at around 14 terms. The effect of including more terms is negligible on the 
results we give in both tables until Nf ^ 8 when perturbation theory can not be regarded as 
reliable anyway. The remaining columns of our tables are the results of this analysis. Clearly 
the effect the 0(l//3o) corrections has is not to move g^^ significantly from the three loop value 
for a large range of Nf. Also for A^ = 14, 15 and 16 the values of f3'{gc^) are not that different 
from the perturbative estimates of f3'{g^'^). 

In conclusion we have produced the leading order corrections to the QCD /3-function in a 



1/A//- expansion which extends the calculation of |12]. Consequently we examined the effect 
they had on the known infrared fixed point in four dimensional /3-function. It transpires that 
perturbation theory is valid for analysing the fixed point when the value of A^ is near the upper 
bound for the existence of gf^ and estimates for a critical exponent were obtained then. It 
would be interesting, though, to compare the values of these exponents with results from other 
techniques such as the lattice, which would be expected to reliably cover the lower part of the 
range. In this case a resummation would be necessary to try and improve the lack of convergence 
which is apparent when the three loop values are compared to the two loop ones. Further, we 
believe it would be useful to repeat our analysis for the supersymmetric extension of QCD in 



large Nf in relation to [23, Once the analogous expression to (10) is available it would be 
possible to study the effect the 0{\/Nf ) corrections have on the infrared fixed point when A'^c is 
large as well as for orthogonal and symplectic gauge groups. As a first step one would need to 
determine which field theory supersymmetric QCD is equivalent to at the d-dimensional fixed 
point and verify, for example, that the correct triple and quartic interactions are obtained in 
the large A^ limit similar to It would be hoped that there is a small set of interactions, as 
in the non-abelian Thirring model, to reduce the amount of calculation that would occur. 
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Nf 


Two loop 


Three loop 


Two loop 
+ 0{l/po) 


Three loop 
+ O(l//3o) 


6 


- 


4.050746 


(0.803122) 


0.517376 


7 


- 


0.782073 


(0.829257) 


0.463261 


8 


- 


0.466021 


(0.858331) 


0.388537 


9 


1.666667 


0.327211 


0.878128 


0.308364 






U.z4ozy / 


U.oOoooz 


U.zooolU 


11 


0.392857 


0.184156 


0.666176 


0.183119 


12 


0.240000 


0.138426 


0.283504 


0.138208 


13 


0.148936 


0.100763 


0.155791 


0.100726 


14 


0.088496 


0.068282 


0.089480 


0.068278 


15 


0.045455 


0.039271 


0.045534 


0.039270 


16 


0.013245 


0.012647 


0.013246 


0.012647 



Table 1. Location of infrared fixed point g^^ for SUc{S). 



Nf 


Two loop 


Three loop 


Two loop 

+ 0{l/Po) 


Three loop 

+ 0(l//3o) 
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81.682972 


(15.663662) 


8.095948 
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5.972522 


(17.792600) 


4.922952 
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2.658882 


(16.790094) 


2.717457 


9 


4.166667 


1.475455 


13.456718 


1.510817 


10 


1.522523 


0.871775 


7.901656 


0.879995 


11 


0.720238 


0.516977 


1.617025 


0.518561 


12 


0.360000 


0.295517 


0.360750 


0.295784 


13 


0.173759 


0.155581 


0.173789 


0.155616 


14 


0.073746 


0.069899 


0.073751 


0.069903 


15 


0.022727 


0.022307 


0.022727 


0.022306 


16 


0.002208 


0.002203 


0.002207 


0.002203 



Table 2. Values of P'{g^^) for SUc{S). 
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Fig. 1. Graphs for anomalous dimension of (Guv) • 
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